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Abstract 

To describe dark energy we introduce a fluid model with no free 
parameter on the microscopic level. The constituents of this fluid 
are massless particles which are a dynamical realisation of the unex- 
tended D = (3 + 1) Galilei algebra. These particles are exotic as they 
live in an enlarged phase space. Their only interaction is with grav- 
ity. A minimal coupling to the gravitational field, satisfying Einstein's 
equivalence principle, leads to a dynamically active gravitational mass 
density of either sign. A two-component model containing matter 
(baryonic and dark) and dark energy leads, through the cosmological 
principle, to Friedmann-like equations. Their solutions show a decel- 
eration phase for the early universe and an acceleration phase for the 
late universe. We predict the Hubble parameter H{z)/Hq and the 
deceleration parameter q{z) and compare them with available exper- 
imental data. We also discuss a reduced model (one component dark 
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sector) and the inclusion of radiation. Our model shows no stationary 
modification of Newton's gravitational potential. 

1 Introduction 

Astrophysical observations (supernovae data [I], [2]) suggest that the uni- 
verse is undergoing an accelerated expansion. This conclusion was drawn 
by interpreting these data in the framework of the cosmological Friedmann 
equations which describe the universe as being homogeneous and isotropic 
on the largest scales (cp P]). Within this framework the origin of the cosmic 
acceleration is attributed to an exotic component, called dark energy, which 
is the source of repulsive gravitation (due to its negative pressure - according 
to the present interpretation). 

But there exist other interpretations of the astrophysical data which do 
not invoke dark energy: 

• Cosmic acceleration could be an apparent effect due to the averaging 
of large scale inhomogeneities in the universe (see [1] and the literature 
quoted therein. For a non-expert explanation see [5]). However, it is 
an open question as to whether this interpretation is in an agreement 
with all available cosmological data (see [6], section 5.3). 

• Modification of the geometric part of the Einstein-Hilbert action by 
replacing the Ricci scalar R by an arbitrary function of it f{R) or 
by introducing higher-order derivative terms (see [7J and the literature 
quoted therein). Some models are based on modified teleparallel gravity 
(see [S]). All such models, however, suffer from having to rely on an 
arbitrary function which cannot be derived from more fundamental 
assumptions. 

Hence we assume that some sort of dark energy is the cause of the cosmic 
acceleration. Before we present our model let us give a very brief critical 
overview of the presently available dark energy models (see also |9]). 
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The simplest model, also called ACDM model, (see any review of dark 
energy, eg [6]) involves the use of a positive cosmological constant A whose 
value has to be determined from experimental data. Its small value (as 
determined by such considerations) causes some problems when we interpret 
A as the energy density of the vacuum (cp JTU], [H])- The most popular 
dynamical dark energy models use instead a scalar field (see the reviews in 
[lOj). However, such models have less predictive power as one can always 
construct a scalar field potential that gives rise to a given cosmic evolution 

Another class of models unifies dark matter and dark energy into a one 
component dark sector. Then the acceleration comes from a new kind of 
interaction within the dark sector. In the case of a Chaplygin gas this inter- 
action is given by an ad hoc assumed equation of state with negative pressure 
(see jllj and the literature mentioned therein). Other models use a complex 
scalar field (see p^] and the literature cited therein) or a phenomenological 
antifriction force which can be understood as a non-minimal coupling of the 
cosmic gas to the curvature |13]. 

In summary; so far we do not have any dark energy model which has been 
derived from fundamental physics [0] . AH known models contain at least one 
new parameter in the microscopic action |14j . 

In this paper we introduce our dark energy model which, on the micro- 
scopic level, contains no new parameters. To do this we start with the well 
known fact that cosmology can be discussed without using general relativity 
as the basic Friedmann equations can be derived from the Newtonian gravity 
(cp. [15j). If we now want to consider some new nonrelativistic particles as 
the cause of cosmic acceleration they must necessarily be massless as the 
usual massive particles always lead to attractive gravitation. The possibility 
of having nonrelativistic massless particles as a dynamical realisation of an 
extended Galilei algebra has already been discussed in some of our recent pa- 
pers ^7\. In the present paper we show that massless particles can exist also 
as a dynamical realisation of an unextended Galilei algebra (a related reali- 
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sation for massless particles moving with infinite velocity has quite recently 
been found by Duval and Horvathy [T6]). 

The existence of nonrelativistic massless particles may appear strange at 
first sight; however, as we show in section 2, such particles possess a modified 
relation between energy and momentum (or velocity) and so they live in an 
enlarged phase space. For this reason we will call these particles 'exotic'. 
Due to the enlarged phase space we have some freedom on how to introduce 
the gravitational coupling for our particles. Here we will do this in a minimal 
way which satisfies the general form of the Einstein equivalence principle but 
which does not use the concept of a rest mass of the gravitating particle. 
This can be stated in the form of the requirement that "a freely falling 
observer does not feel any effect of gravitation" [T8]. This minimal coupling 
has the important property that, in a many exotic particle system, it leads 
to a dynamically generated active gravitational mass density of either sign 
which can then be a source of the gravitational field. Then a fluid mechanical 
generalisation of this model can serve as a new model for dark energy. 

A further extension, to a two-fluid model, including baryonic- and dark 
matter besides dark energy, then leads, using the cosmological principle, to 
Friedmann-like equations for the cosmological scale factor a{t). The solutions 
show a deceleration phase for the early universe and an acceleration phase 
for the late universe. 

Furthermore, we show that our model allows for a one component descrip- 
tion of the dark sector on large {ie cosmological) scales. The choice between 
these two possibilities has to be made by comparison with the observational 
data on galactic scales. 

We have also looked at the influence of our dark energy sector on local 
systems. We show that, in particular, it does not lead to the modification of 
Newton's gravitational potential. 

The paper is organised as follows. In section 2 we present our nonrela- 
tivistic massless particle model coupled minimally to gravitation. In section 
3 this model is generalised further and extended to a two component fluid 



4 



model for matter (baryonic and dark one) and dark energy. In section 4 
we show that we do not run into instabihties with our modeL In section 5 
we describe some solutions, which satisfy the cosmological principle, of the 
corresponding fluid dynamical equations. These solutions show a decelerat- 
ing universe at early times and an accelerating one at late times. We will 
see that on the cosmological level our model is an effective one component 
model for the dark sector. In section 6 we include radiation. Observational 
consequences discussed in section 7 are: 

• The prediction of the Hubble parameter H{z)/Hq and of the deceler- 
ation parameter q{z), having fixed two integration constants by using 
measured cosmological parameters. 

• The proof that Newton's gravitational potential requires no stationary 
modification. 

Some technical details are given in appendix A. In appendix B we spec- 
ulate on a relativistic generalisation of our nonrelativistic particles describes 
tachyons. We conclude with some final remarks (section 8). 

2 Nonrelativistic massless particles and their 
gravitational interaction 

In our second paper in [T7] we have introduced the Lagrangian 



where, in the three dimensional {Ui) are the components of spatial 

position (velocity) of a particle and Pi (g^) are the components of the cor- 
responding momenta. We use Euclidean metric and Einstein's summation 
convention with z = 1, 2, 3. An overdot represents a time derivative. 

The Lagrangian ([1]) leads to a dynamical realisation of the acceleration- 
extended Galilei group in any dimension with one central charge (k) for a 




(1) 
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non-interactive massless particle. Without the last term in ([T]) we have a 
dynamical realisation of the Galilei group without any central charge {ie 
without any free parameter). 

To show that we note that when k = oo, the equations of motion that 
follow from ([T]) are 

ii = Vi, Pi = 0, qi = -pi, iji = 0. (2) 

These equations correspond to the canonical Poisson brackets (PBs) 

{Xi, pj} = Sij, {yi, qj} = Sij, (3) 

which can be derived from the Hamiltonian 

H = PiVi. (4) 

If we now introduce the conserved Galilean boost generator Ki which is 
given by 

Ki = Pit + Qi (5) 

we find that 

{P^, K,} = 0, (6) 

which clearly shows that we are dealing with a massless particle. 

Going the other way round, ie by starting with m = 0, with ([6]) as a 
requirement, it can be shown that the Lagrangian 

L = Pi{xi - Vi) + Qiiji, (7) 

defined in a 12- dimensional phase space, is the minimal one [20j . 

Furthermore, we note that the conserved angular momentum is given by 

J = xxp + yxq (8) 

and that the Poisson brackets of p, K, J and H build the unextended Galilei 
algebra. The presence of the second term in ([8]) shows that our particles 
possess a nontrivial spin. 
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To couple this particle to gravity we start with the general form of Ein- 
stein's equivalence principle. In a nonrelativistic context this can be stated as 
follows: locally, ie at each fixed space point x, a gravitational force —\/(j){x, t) 
is equivalent to a time-dependent acceleration h{t). The only known equa- 
tion of motion for the particle trajectory x{t) satisfying this form of the 
equivalence principle is given by the Newton law: 

Xi{t) = -di(j>{x{t),t) (9) 

because is invariant with respect to arbitrary time- dependent translations 

(cp m) 

Xi ^ x[ = Xi + ai(t) (10) 

if 0(a;, t) transforms to 

(j)'ix',t) = (f){x,t) - di{t)Xi + h{t). (11) 

Hence considering </)(x, t) as an external gravitational field we can take 
its interaction term with our particle Li^t in the form: 

Lint = qidi(f){x{t),t) (12) 

Clearly, with this term, the equation of motion for Xi is given by iQ and 
the second equation in ([2]) becomes pi = qkdkdicj). Then our system is in- 
variant under arbitrary time- dependent translations ffTOj) where transforms 
according to f|TT]) with qi and Pi being invariant. 

3 Two-fluid dynamics 

In this section we consider a two-fluid cosmological model where one fluid 
component M consists of massive matter (baryonic and dark one) and the 
other fluid D consists of the exotic massless particles, introduced in the pre- 
vious section and representing dark energy. The only interaction considered 
within the fluids and between them is gravitational. 
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3.1 Lagrange picture 

First we generalise the dark energy model introduced in the previous section 
to the continuum case by introducing comoving coordinates ^ E [21], add 
continuous massive matter with its standard gravitational interaction and 
use the usual Lagrangian for the gravitational field. 
Then our Lagrangian is given by 

L = Lm + Ld + L^, (13) 

where 

LM = mj £i (yf (xf - \yf) - </>(f^^t)) , (14) 
where m is a mass parameter giving f|T^ the correct dimension, 

Ld = j d'i (p,(i;f -yf) + gfyf + g.5,0(f^,t)) (15) 



and 



StiG 

In these expressions all phase space variables are functions of ^ and t, ie 
^ x^^i^^t) etc. 

Note that both Ld and Lm are invariant, up to a total time derivative, 
under the transformations fllOtilip . 

The equations of motion (EOM) corresponding to L are given by 



M sector 



if = 



yf = -d,<j>{x'\t) (17) 



D sector 



if = 



Qi = -Pf (18) 

yf = -9,0(x^,t) 
pf = qkdkdi(j){x^,t) 
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• sector 

A0(x,t) = 47rG j d^^(mS{x-x^{i,t)) + q,{tt) di5{x - x^ {(,1))^. 

(19) 

The last term in fll9p represents a dynamically generated active gravitational 
mass density. 

3.2 Eulerian picture 

In the Eulerian picture the dynamics of the fluid is described in terms of 
X and t dependent fields: particle number density n{x,t), velocity Ui{x,t), 
momentum pi{x,t) and pseudo- momentum qi{x,t). 

Assuming uniform distribution in ^ the Lagrangian phase space variables 
are transformed to the Eulerian fields by 

n{x,t) = j d^^5%x- x{i,t)) (20) 

and 

n(f,t)p,(f,t) = j d'^piii,t)5%x-xii,t)) (21) 

and an analogous expression for Ui{x,t) (in the expression above replace 
Pi{x, t) by Ui{x, t) and Pi{C,, t) by yi{C,, t)). Similarily for qi{x, t). In fact, fl2Tl) 
holds for any function of relevant variables. 

To derive the EOM in the Eulerian picture we follow the standard proce- 
dure (cp. [21]) and obtain from fll7til9p by using (I2UII21I) the corresponding 
equations in the Eulerian picture: 

dtn^{x,t) + dk{n^ut){x,t) = 0, (22) 

where A = (M, D), ie the continuity equations for the particle number den- 
sities ra^'^ and n^, and from (|T9l) the Poisson equation for the gravitational 
field 

A0(f,t) = 4nG (p^ + a,(n^g,)) , (23) 
where the mass density p*^ is defined by p^^ := mn^ . 
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Note that the last term in (1231) represents the dynamically generated 
active gravitational mass density of the dark-energy fluid. 
We have, in addition, the following Euler equations: 

«f = -84 (24) 

(from the second equation in (fT7|) ) and from the third equation in (fTSl) 

Dfu^ = -04, 

where we have defined Df = dt + ufdi. 

Suppose now that -uf^ and obey the same initial conditions. Then 
(!24|) shows that uf = uf^ = Ui ie (|24|) becomes one universal Euler equation 
valid for all fluid components. 

DtUi = -84- (25) 

Finally, the second and fourth equations in (fTSjl give 

DtQi = -Pi, DtPi = qkdidkct). (26) 

Looking at ( 125)1261) we note that, in contrast to standard fluid mechanics, the 
two vector fields p(a;, t) and ^(x, t) are not parallel to each other. 

3.3 Symmetries 

First we note that all our EOM (12211261) are obviously rotationally symmetric. 

To consider other symmetries we observe that if we perform an infinites- 
imal time dependent translation 5xi = ai{t) we see that 

Sui{x,t) = di(t) - ak(t) dkUi{x,t) 

S(f){x,t) = -di{t)x, + h{t) - ak{t)dk(j){x,t) (27) 

and 

5C{x,t) = -akit)dkC{x,t) 
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where C ^ i^^yPi^Qi)- Thus the EOM are invariant under such translations 
and so, locally, the general form of Einstein's principle of equivalence is sat- 
isfied as in General Relativity. Moreover, as shown recently by one of us 
(PCS), we obtain, when neglecting the massive matter part, as symmetry al- 
gebra, the expansion-less conformal Galilei algebra with dynamical exponent 

z = im- 

3.4 Stress tensor and pressure 

To see that our massless particle fluid may, indeed, represent dark energy we 
show now that the pressure can be negative. To do this we consider the local 
momentum conservation 

dtPi{x,t) + djT,,{x,t) = 0, (28) 

where, in our case, the momentum density Pi and stress tensor Tij are given 
by 

Pi{x,t) = {npi){x,t), T,j{x,t) = iP,Uj)ix,t) + V5,j (29) 

and V is the pressure. 

Note that the stress tensor is not symmetric. The deeper reason for that 
is the presence of a spin part in the conserved angular momentum (see ([8]) 
and cp. [22]). 

To have the system as simple as possible we consider first a one- dimensional 
self-gravitating massless particle fluid. The EOM for the momentum field 
p{x,t) then, due to the Poisson equation (fT9l) . becomes 

Dtp = A7iGqd,{nq). (30) 

Using the continuity equation (!22l) . together with (!30l) we obtain for (p8|) 

dtP + {Pu - 2nG{nq^)) = (31) 

ie we get the negative pressure 

V = -2710 {nqf. (32) 
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For a three dimensional case we obtain a result similar to 0321) if we assume 
that the vector field ng is irrotational (see section 5), ie 

{nqk){x,t) = dkh{x,t). (33) 

Then, from (fT9|) . we get 

di^ = AnGnqi, (34) 
leading to the EOM for the momentum field Pi 

DtPi = AnGqkdiinqk) (35) 

and finally, analogously to (!3T|) . to the pressure 

V = -2TiG{nqkf. (36) 

If nq is not irrotational then the stress tensor Tij is more complicated than 
that given in fl29|) . 

3.5 Nature of the gravitational mass 

According to the Poisson equation ( l23l) the active graviational mass density 
of the dark energy fiuid is given by 

= diinqi). (37) 

This result should be compared with the expression, from general relativity, 
for the perfect fiuid 



= [P" + f ) . (38) 

where c^p^ is the energy density in the rest frame of the fiuid and V - 
the corresponding pressure. Neglecting the massive matter component and 
considering the particular case of an irrotational field nq (see ( 133|) ) we find 
from the Lagrangian (11511161) and the transformations (120112 ip that the energy 
density S^{x,t) is given by 

= npiUi - 2nG{nqif (39) 
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and so 



= -'J^^^. (40) 



Then, with V given by fl5Bl) . we obtain for 



= (41) 



The reason for the difference between pil) and fl37j) is obvious: in our 
nonrelativistic massless particle model there is no place for the velocity of 
light c! Furthermore, we note that fl371) is not built from the components 
of the energy-stress tensor and, in fact, fl37|) is the simplest Galilei invariant 
expression in our model with the dimension of mass density. 



4 Stability considerations 

The Hamiltonian corresponding to the Lagrangian f|T5|) of the dark energy 
fluid component is linear in the momenta and thus not bounded from 
below. This property is well known for any higher-derivative Lagrangian 
(theorem of Ostrogradsky [23]) and it arises in our case because our one- 
particle Lagrangian considered in section 2 may be understood as the limiting 
case of the higher-order Lagrangian ([T]), which in the configuration space 
takes the form L = |xf (see [I7j2). However, for a free particle, this does 
not concern us too much, as the Hamiltonian (jl]) can always be transformed 
to a positive quadratic form by a complex-valued canonical transformatiorl^: 

Pi = iap[ - bq-, yi = -hq[ - iap[, 

where a, b are arbitrary real numbers, leading to 

H = a'pf + iPqf. (43) 



^Similar ideas can be used to demonstrate the absence of ghosts in the Pais-Uhlcnbeck 
model. This was recently shown, in a different way, by Bender and Mannheim [24j 
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Thus all serious problems (like, eg the collapse in the classical case or a 
nonunitary time-development in the quantum case) may arise only in the 
presence of interactions (cp [SB]). 

Another possible instability of interacting theories containing negative 
energy involves the spontaneous decay of any state into a collection of positive 
and negative energy particles (cp. |25j)- This instability is excluded in our 
model due to the particle number conservation fl22|) . 

To find out what happens in our fluid model it is sufficient to consider 
only a self-gravitating massless particle system, ie to neglect its dark matter 
component. To treat things as simple as possible, in the following we consider 
only a one-dimensional system. In subsection 4.1 we show that a two-particle 
system at zero energy can, indeed, collapse. When, in the next subsection, we 
generalise these considerations to the continuum, ie a fluid dynamical case, 
we find that in this case the collapse does not take place. This reassures us 
in our belief that the three dimensional case is also collapse free. 

4.1 Two-particle case 

Specialising the EOM (ITSl) . for = 1, to the two-particle case we obtain for 
the relative motion: 

x{t) = A7rGq{t)S{x{t)) (44) 

and 

q{t) = 2nGq\t)6'{x{t)) (45) 

where we have defined the relative variables (the indices 1,2 label the two 
particles) 

X := xi - X2, q ■= qi - q2- 

The variables of the two-particle centre, R := ^{xi + X2) and Q := \{qi + q2) 
satisfy 

Q = 0, and R = -A7rGQS{x). (46) 
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To obtain these equations we have used the fact that, for a generic x & R^, 
dx4>{x,t) is given, due to (fT9l) . by 

2 

d,(P{x,t) = AnGj2 Six - xa{t))qa{t); (47) 

a=l 

- here we have taken the particular solution Q{t) = of fj46|) . 
Then the energy of the two-particle system is given by 

E = + nGq^S{x). (48) 

A solution of the EOM fj44l) . fH5|) . for vanishing energy E, is clearly given 

by 

X = 27rGXq\ q = X'^ix), (49) 

where A is arbitrary and needed for dimensional reasons. 
If we now take A > so that 

x{t) > 0, for all t e (50) 

and choose 

x{t) < 0, x{t) > for t <to and x(to) = (51) 
we obtain from (H^ that 

q{t) = q,{l-e{t-to)) (52) 

and 

(2nG\q^ for t < to 
x{t) = \ (53) 

[O for t > to 

ie the two particles collide at t = to aiid stay together for all later times. 
This is a collapse situation. 
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4.2 Hydrodynamic case 

In the Eulerian picture for d = 1 the analogue of (H7|) is now 

d^(f) = AixGnq (54) 

and so the hydrodynamic EOM obtained from fl22p . fl25l) and fl26p take the 
form 

dfU + dx{un) = 0, (55) 
Dtu = -AnGnq (56) 

and 

D^q = -AirGqd^inq). (57) 
In analogy to ( l39l) the energy density S{x,t) is given by 

^ = n (-(Ag)M - 2TxGnq^) . (58) 

Next we proceed as in the two particle case. We make the ansatz 

u = 2nG\q^. (59) 

Then from (l56l) we find that 

Dtq = -X-^n (60) 



which, together with fl59|) . demonstrates the vanishing of the energy density 
S. 

Then, as can be easily checked, the unique solution of the remaining EOM 
dSSD, ([57D and ([59]) is given by 



n{x,t) = no (= const), q{x,t) = — (^■'■) 

A 



demonstrating a collapse-free situation. 
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5 Cosmological solutions of fluid dynamics equa- 
tions 



In order for the universe to be homogeneous and isotropic on large scales we 
require, as usual, that 

= n^it) (62) 

and 

Ui = ^Xi, (63) 
a[t) 

where a{t) is the cosmic scale factor. 
Then ((251) tells us that 

dicj) = Xiipit) (64) 

with 

, , a 
^{t) = —. 

a 

Putting ( l63l) and (Ell) into the second equation in ( 126|) gives us 

a 

DtPi = -Qi-- (65) 
a 



To solve the first equation in ( 12611 and ( 16511 we make an ansatz 

Qi = fq{t)xi, and Pi = fp{t)xi. (66) 
Then, using (I63ti65p . we eliminate fp and get 

/, + 2-/, = 0, (67) 

(X 

which can be integrated once giving us 

fq{t) = -J^ with 13 = const. (68) 

Furthermore, with (IB^ and (IB^ the continuity equations (12^ can be inte- 
grated as usual giving us for p*^ and 

M M 
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and 

where M and D are positive constants. 

Inserting flM|) and fl69|) into the Poisson equation fl23|) we get from (|66|1 

-a = ^(M + 3Z}/,), (71) 

where fq should be taken as a solution of (l68l) . Eq. (ITTjl is one of our 
Friedmann-like equations. 

We should now distinguish two cases: 

• (3 = which implies /q=const. Then eq. (17T|l gives us that 

M 

a > for any t if fg< -— (72) 

ie we obtain an accelerated expansion for all times (this contradicts the 
known cosmological facts). 

• /? 7^ 0. Then putting (168!) into (iTTj) we get 



P 

Integrating once we find 



(M + 3D/,). (73) 



a = ^iM+^Df,)+c, (74) 



Multiplying (1741) by fg and using (!68|) on the l.h.s we obtain 

^ ^^(M+^D/,) + ci/, (75) 



a2 /? 
which after integration gives us 



^ = ^ + ^/.) + cifg + Co, (76) 
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where Cq and Ci are integration constants. 

Let us now discuss, given (J76il . the behaviour of fq as a function of the 
scale factor a. 

Performing the transformation: 

/. - 9{a) ■■= h + ^ (77) 
we arrive at (redefining cq and Ci) 

g\a) + cig{a) + cq (l - -) = 0, (78) 
where we have defined 

Let us look now at the solution of ( ITHI) with the constants cq and c\ being 
positive, co,i > 0. First of all we note that the scale factor a may serve as a 
measure of time due to d > (expanding universe). 

For a < a* we have gia) > and so, due to ( 17T1) d < 0. So, for a < at, 
we are in the deceleration phase of the early universe. On the other hand, 
clearly, for a > we have g{a) < and then, due to (1711) . d > 0. So, for 
a > at we are in the acceleration phase of the late universe and we see that 
at defines the transitional scale factor at which the deceleration stops and 
the acceleration takes over. It can easily be seen that the condition co,i > 
is also necessary to obtain these results. 

Next we observe that by differentiating (ITHj) with respect to a we have 

g'ia) = -^r^T^' (SO) 

a"^(ci + 6g'^[a)) 

where ' denotes the derivative with respect to a. If we now put f lHOj) into fl68l) 
we get 

a = Jj£l±im (81) 

Co at 

thus showing that, for d > 0, we need (3 < 0. 

Eq. (1811) is our second Friedmann-like equation. Note that the first 
Friedman-like equation flTTj) is a consequence of the second one ( IHTj) if g{a) 
is a solution of the cubic equation fl78l) . 
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To integrate fl8ip we need the explicit form of g{a). To obtain g{a) we 
note that g{a) is the real valued solution of the cubic equation (ITHl) . This 
solution is given by 

g{a) = u+{a) + u_(a) (82) 



with 



where 



u±{a) = ( -| ± 



Ci\3 /g\2 



(83) 



(- 

Then, from flHTl) we find that 



t-to = ^ / da ; (84) 
' ci + 35(^(a) 



with (7(a) given by ( l82l) . 

In the Appendix A we present a detailed discussion of the evaluation of 
fl8^ in terms of roots of flTSl) . As our final results are not very transparent 
let us mention here some asymptotic results: 

• At large a, ie a ^ at, a grows linearly with t. This follows from the 
observation that at large a, q goes to Cq and so the integrand in (IMI) 
becomes independent of a. 

• At a very close to at we get from flTSl) that 

g{a) ~ — —{a -at). (85) 

Ci at 

Then, by choosing to as the time at which a = aj we obtain from (IMll 

, , -i , l + i-f{a{t) - at) , . 

t-to ^ — log — -, 86 

d 1 — i'y[a[t) — at) 

where 7 := #a and S := 

c( at cf a^ 

Inverting (186!) and taking the first terms of the power series expansion 
in t — to we obtain 

a{t) ~at = ^(t - to) + ^(t - to)' + 0{{t - to)'). (87) 
Coat Coat 
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Considering (15^ with (ITTi) and (ISO]) it is easy to see that higher order 
corrections to (l85l) do not change the first two terms in the expansion 

• For small a ie for a <^ we obtain from (1781) 

- (^fi)* (88) 

leading to, due to (18^ with to = and a(0) = 0, 

a(t) ~tt (89) 

thus showing that the combined effect of matter and dark energy at the 
early times differs from the behaviour of the matter dominated universe 

2 

for which a{t) ~ ts. 

Note that this result (l89l) is exactly the scale invariant solution for a{t) 
corresponding to the dynamical exponent z = 5/3 [20] . 

5.1 Dark sector with one or two components? 

In section 3 we made the usual assumption that the dark sector possesses a 
two-component structure. However, our results for the Friedmann-like equa- 
tions (ITTil . (1ST]) and for the equation (1751) determining g{a) are all independent 
of the constant M, defined by (1^ . Thus, as long as we do not compare (ITTi) 
and (18 ip with the original Friedmann equations which would be physically 
senseless due to the different nature of the gravitational mass in our model 
(see section 3.5), it is sufficient to keep the dark energy fluid (now to be called 
"dark fluid" ) as the only component within the dark sector. In our case, this 
dark fluid takes over the role of dark matter and dark energy, at least on 
large scales, like in the cases of the Chaplygin gas or the complex scalar 
field [T2j. The baryonic component, which corresponds to about 4% of the 
energy of the universe, is negligible on these scales. 

To describe the universe correctly, at the scale of galaxies, the dark fluid 
must behave like dark matter; ie exhibit attractive gravitation at local scales 
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(see [12] and the literature cited therein). This point still has to be examined 
in more detail. 



6 Cosmology including radiation 

Including radiation (photons), and also massless neutrinos, within our frame- 
work, would require a full relativistic treatment. However, what we really 
need here is somewhat less ambitious. For the cosmology as outlined above 
we need a description of radiation as a nonrelativistic fluic|^ component R 
with an equation of state parameter (defined as the ratio of pressure and 
energy density) [6] 

co^ = i. (90) 

To get the required result we follow McCrea [27] and Harrison [28] who 
extended Newtonian cosmology by taking pressure into account. For a homo- 
geneous and isotropic universe we have therefore to add to our hydrodynamic 
equations the continuity equation for the radiation energy density c^p^ 

pR + 4- = 0, (91) 

Qj 

whose solution is given by 

where i? is a positive constant. Furthermore we must change the Poisson 
equation (1231) by adding to its right hand side the active gravitational radi- 
ation mass density 2p^(t) leading to 

A0 = 47rG(p*^ + diin^Qi) + 2p^). (93) 

From fl93l) we conclude that the first Friedmann-like equation flTTl) now 
becomes 

G 2 /? 

-d= -{M + 3Df,+ — ). (94) 



^Note that within a hydrodynamic description of radiation the velocity field at a point 
X is an average over all directions of radiation velocities whose modulus is therefore less 
than c (it might even be small when compared to c). 
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Unfortunately, when i? 7^ 0, it is not possible to integrate analytically the 
coupled system of differential equations and (|68ll . Nevertheless, we can 
conclude, as usual, that at very early times the last term in fl94l) dominates, 
ie the universe is radiation dominated. In the following, we will consider, as 
we have already done in section 5, the universe only for the later times, ie 
when the last term in flMl) is negligible. 



7 Observational consequences 

Our exotic massless particles possess no non-gravitational interaction, nei- 
ther with the particles of the Standard Model nor with the dark matter 
particles. Thus their existence can only lead to observational consequences 
at cosmological scales (see sections 7.1-2) and, perhaps, also at local scales 
(see section 7.3). 

7.1 Predicting the Hubble and the deceleration pa- 
rameters from our model 

To calculate the Hubble parameter H in our model we introduce the redshift 

z by 



and then consider 



a = ^— (95) 

l + z ^ ^ 



H{z):=-^{z). (96) 



Then from (17T!) we find that H{z) is given by 



Hiz) = ^ fl + 3^ 



1 

"l + z 



(97) 



acoOt \ Cl 
Next we define Hq as Hq := H{z = 0) and so find that 

H(z) 1 + 3^^i±^ 

Cl 
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In a similar way we see that the deceleration parameter q{z) defined as 

a 



q{z) 



aH^z) 

is given by 

f^n , n(n\ / n'^inW ^ 

(99) 



Qat Co g{a) ( ^9^{a}^ ^ 

q{z) = 1 + 3 

a 2 „2 \ Ci 



1 + 2 



Cf Ci 

Note that both are functions of only at and oi k = Hence to determine 

"-1 

them we need two experimental data. 

Before we try to determine q{z) and H{z) let us observe that there are 
a few things we can say about their behaviour for any values of the two 
parameters. First of all, we easily see from (188|) that at large z, ie for z 3> 



silt.) 



and so from (!99|) we get that, for z'^ Zt^ 



1 • ■ 1 



3 



Co 



3 



(100) 



q{z) ~ I (101) 

Moreover, /i(z) is monotonically increasing. To see this we take (p7|l and 
note that 

h{z) = k{l + 3f{z)) (l + z), (102) 

where we have defined 



g{z) := Ci 2 g{ 



1 



■l + z' 

and, similarily, the overall positive constant k. Then fISUl) is equivalent to 



Then 



In addition, from fllOOl) and fll02p we see that, for z ^ zt-, 

h{z) ~ ^t. (105) 
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7.2 Estimation of H{z) and of q{z) 

To obtain our 'predictions' for h{z) and q{z) we use the data from the first 
reference in [29j. They give us g(0) = —0.57 and Zt = 0.71, both with small 
errors which we do not mention here as the curves we will show here depend 
very little on the exact values of these parameters. These values are obtained 
by fitting the matter part of the ACDM model to observational data. 
g(0) then serves to determine the constant k, in our model. We will give the 
curves obtained with these values subscript S. 

We can also use the model independent values from the other two ref- 
erences in |29]. The data from the paper by Cunha are g(0) = —0.73 and 
Zf = 0.49 and from the paper by Lu et al g(0) = —0.788 and Zt = 0.632. The 
curves corresponding to them will carry the indices C and L respectively. 
Note that the values of k for the three cases are k, = 0.8667 (S), 0.977 (C) 
and 1.166 (L). 

In fig. 1 and 2 we plot our predictions for H{z) and q{z) respectively. We 
have normalised H{z) to its value at 2; = 0, so in fact, our plots are of h{z) 
In fig. 3 we present the corresponding values for g{z). We note that all 3 
cases are quite similar. 

We have also attempted to compare our results to the experimental data 
given in table 1 in 130] • The results given there have large experimental errors 
and are given only for a few values of z. Hence they will not be too conclusive 
or reliable. However, to perform any comparison we need the value of Hq. 
We can, of course, take this value from [2]2. There we find Hq = 70. 5. Hence 
in fig 4 we present our data (with the normalisation fixed by Hq = 70.5) 
and compare them with the experimental data (constructed from the data 
in [30]) corresponding to the experimental data + 1 standard deviation error 
(called 'maximum') and -1 standard deviation error (called 'minimum'). We 
note the general agreement and so we are heartened by this result. 
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Hubble function h(z) 



our prediction_S 
our prediction_C 
our prediction_L 



Figure 1: Our prediction for h{z) 




Figure 2: Our prediction for q{z) 
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Figure 3: Our prediction for g{^z) 



Comparison with experimental data 
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Figure 4: Comparison with experimental data;i7 (^) in units km sec ^ Mpc 
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7.3 Influence on local systems 

Here we look at the problem of how a two-body system, bound by the stan- 
dard Newtonian potential, may be affected by the dark sector proposed in 
this paper. To study this we consider two different mechanisms: 

• The effect of the dark fluid at cosmological scales giving rise to an 
additional time- dependent term for the two body potential 

50(r,t) = (106) 

The equations for the two-body relative motion then takes the form: 

r = -r — -^r, (lU7j 
a r-^ 

where /i is the reduced mass. 

As we do not have the explicit form of the time dependence of the scale 
factor a (if:) we use instead a as a measure of time. Then ( 1107P leads to 
the following differential equation for f{a): 

r d + r d = —r — -^r (108) 

a 

or using the Friedman-like equations (!7T!) and (IHTI) . we obtain 



[ci + 3^2^a))2 ?,DG 



r {a)(3' ' " —g{a) r (a) = (109) 



-g{a)r{a) ^--r(a). 



where g{a) is given by (182!) and a prime denotes differentiation with 
respect to a. 

To solve (I109p numerically we would have to know besides the con- 

3 

stants at and cq/c^ known from 7.1, also the values of constants (3 and 
Cq. Recent estimates of the effects caused by 6(f){r,t) in the case of a 
constant < — 1 [31] have found observable effects on a time-scale 
given by billions of yearj^. We expect similar results for our model. 



■^For consideration of more general astronomical structures see [32] and the literature 
cited therein. 
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• The other issue involves a possible modification of Newton's gravita- 
tional potential by a local, stationary dark energy fluid. To study this 
we consider a point mass m located at x = 0. We will show that the 
corresponding stationary dark energy flow leads to a vanishing extra 
gravitational mass density 5j(p^gf ) and so there is no extra contribu- 
tion to 0(r). To see this we consider the D sector of our equations of 
motion given in subsection 3.2 for the stationary case. They become 

dkin'^uu) = 0, (110) 

UkdkUi = -di(j) (111) 

Ukdkqf = -P? (112) 

Ukdkpf = q^dkdict> (113) 
together with the Poisson equation 

A(/. = 47rG(m5(f) + d,{n^qf)). (114) 

Then we use f lll2p and fillip to eliminate and dicp in f lll3p and 
obtain 

Uk dk ui di qf = dk ui di Ui. (115) 

Looking at flllSp we note that it implies that qf has to be proportional 

to Ui 

q? ~ Ui (116) 
and so, due to fill OP . we obtained the desired result ie 

dk {n'^q^) = 0. (117) 

8 Final Remarks 

Given that there are already many dark energy models, what are the reasons 
why we have introduced a further one? The reasons are twofold: 
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• There are no free parameters in the microscopic formulation of our 
model. 

• Our model introduces new physical ideas in the form of nonrelativis- 
tic massless particles whose minimal coupling to gravity leads to the 
generation of an active gravitational mass density of either sign. 

This last point poses the question about the relation of these new physi- 
cal ideas to Newton's and Einstein's theory of gravity. As our particles are a 
dynamical realisation of the unextended Galilei algebra they fit into the gen- 
eral scheme of nonrelativistic physics. The gravitational coupling, satisfying 
Einstein's equivalence principle, leads to the same equation of motion (8) in 
configuration space as in the massive case. Thus we can consider our model 
of a gravitationally coupled, nonrelativistic massless particle as an extension 
of Newton's theory of gravity. 

However, it seems to be not possible to obtain our model as a nonrel- 
ativistic limit of a relativistic model. Massless relativistic particle models 
possess conformal Poincare symmetry leading, in the nonrelativistic limit, 
to conformal symmetry ^7\i, ie z = 1. In our case, for a one-component 
dark sector, = §• In Appendix B we have speculated that the relativistic 
generalisation of our Galilean massless particles are tachyons. However, it 
may be that we are here in a situation similar to Hofava gravity [33j; ie we 
have nonrelativistic symmetry in the ultraviolet limit (small t) and approach 
General Relativity only in the infared (large t) limit. However, to have such 
a picture we may have to modify our model. This is a challenge for further 
research. 

As a drawback of our model one can consider the existence of additional 
dimensions in phase space. However, such a case is already well known 
from the related case of nonrelativistic massless fields (Galilean electromag- 
netism) in which the lagrangian formulation requires the introduction of aux- 
iliary fields [34j. In our case the additional degrees of freedom lead in the 
Friedmann-like equations to undetermined constants which are integration 
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constants along the additional phase space dimensions. The question then 
arises as to whether these constants can be determined a priori by some 
physical arguments. This point is currently under investigation. 



9 Appendix A 



Here we demonstrate that the integral (184|) can be calculated in a closed 
form. 

First we note that due to fISUD we have 



:ii8) 



Coat / da(ci + 3g^{a)) = — da a g'{a). 



Next we change the integration variable a g{a) and use (I75|) to rewrite 
the right hand side of (11181) as 



{coatf j dg {g^ + Cig + Cq) 



We define the roots of the cubic equation 



as gi. They are given by 
gi = v+ + g2 

with 



g + cig + Co = 



7^ 1 7^ w3, gs = g^, 



Co 



Next we perform the decomposition 

3 



Cl 



3J ^ V2 



Co 



(g^ + cig + Co) ^ = Ylig -gi) ^ = ^^^^ ^ ^' 

i=l i=l 

where Oj are given by 

o-i = i{9i- 9i+i)i9i- 9i-i)y^ ■ 



;ii9) 



(120) 



;i21) 



;i22) 



(123) 



(124) 
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Here i = 1, 2, 3 and cyclic permutation is assumed. 

Putting all this together we perform the integration in flllQp and obtain 

/I ^1 
ci + 3g^{a) g{a) - g, 

ttittj g{a) - gi 



2 Coat V ^log^^ ■ 

gi- 9 J g{a) - gj 



Clearly oi = a* and 03 = 03. 

10 Appendix B 

Here we discuss a possible relativistic correspondence of the nonrelativistic 
massless particles introduced in section 2. Clearly, they cannot correspond 
to either massive or massless relativistic particles. However, they could cor- 
respond to tachyons which can be seen as follows: 

• The relativistic generalisation of the equations of motion ([2]) are given 
by the derivatives of the corresponding four-vectors with respect to the 
relativistic parameter r: 

X,, = Vix, Pix = 0, = -p^, y^, = 0. (126) 

• From the second and fourth equations in 01261) we see that 

p^y^' = const. (127) 

However, in order to reproduce, in the non- relativistic limit, the energy 
relation (jl]) the constant appearing on the right hand side of fll27p must 
vanish, ie we must have 

P.y'' = 0. (128) 

• From 01281) we see that 



P.P" = - f =<-ll - -]f < 0. (129) 
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